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Experiments at high luminosity electron-positron colliders allow to study the kaon and the pion form factors 
in the time-like region up to high energies. Also the kaon and the pion pair production at and around the narrow 
resonances J ftp and '0(25') can be investigated, with the interference between electromagnetic and hadronic 
amplitudes as one of the most interesting phenomenas. Parameterisations of charged and neutral kaon as well as 
pion form factors, which lead to an improved description of the data in the region of large invariant masses of the 
meson pair, are presented. They are implemented into the Monte Carlo generator PHOKHARA, together with 
the hadronic couplings of charged and neutral kaons to J/t/i and '0(25). The physics case as well as details of the 
implementation are discussed. 



1. Introduction 

During last years meson factories, using the ra- 
diative return method |3I4) . contributed a lot to 
the error reduction of the hadronic contributions 
to the muon anomalous magnetic moment and 
the running electromagnetic fine structure con- 
stant [iPj . 

An important tool in these analyses is a Monte 
Carlo generator which simulates all the measured 
reactions. To meet these requirements, the gen- 
erator EVA was developed |4l5j . which is based 
on leading order matrix elements combined with 
structure function methods for an improved treat- 
ment of initial state radiation. To improve the ac- 
curacy of these simulations the complete next-to- 
leading order (NLO) QED corrections were eval- 
uated [617] and implemented into the generator 
PHOKHARA |8|9|10|11|12|13|14|15|16|17] . For a 
recent review of theoretical and experimental re- 
sults see e.g. [18]. 
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In the last paper, resulting in the implementa- 
tion in the PHOKHARA generator [17 , the pos- 
sibility of studies of the narrow resonances J/ip 
and ^'(25') at B-meson factories was investigated 
in details profiting from the previous work on 
the narrow resonances [19j . In this contribution 
the results of |19I17| are summarised and com- 
plemented with additional material not presented 
there. 

The new models of pion and kaons form fac- 
tors are presented in Section [21 In Section [3] 
the modelling of the narrow resonances is shortly 
sketched. In Section 2] the implementation of the 
narrow resonances in PHOKHARA is described 
and a discussion of the importance of the FSR 
NLO radiative corrections is presented. A short 
summary follows in Section [5j 



2. Modeling of the pion and kaon form fac- 
tors 



For 



a 



realistic 



generation of the processes 
e+e^ TT+TT^ + photons, e+e^ — K^K~ + 
photons and e+e^ — > K'^K^ + photons mod- 
els for the electromagnetic pion and kaons form 
factors are required. In PHOKHARA 5.0 |14I13| 
and 6.0 [TS] models presented in [20 were imple- 
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merited. They were published before the CLEO- 
c measurement of the form factor in the vicin- 
ity of the '0(25') resonance [i2I and underesti- 
mate the experimental results significantly. Sim- 
ilarly the model predictions at J/ij) are lower 
than the pion form factor calculated in [22] from 
B{J/ip — )■ 7r+7r~) and B{J/ip — e'^e~) decay 
rates. In [T^ an effort was made to get better 
description of the form form factors, especially at 
high meson pairs invariant masses. In combina- 
tion with the J/ip and ipi'2S) hadronic amplitudes 
parameterisations, obtained in |19| on the basis of 
new experimental data, it allows for a description 
of the pion and the kaon pairs production reflect- 
ing all existing experimental data. 

The models proposed in [17] are generalisations 
of the models used in [20] and for both, pions 
and kaons, they are inspired by the dual QCD 
model of the pion form factor [23]. The infinite 
towers of p, w and (j) radial excitations, present 
in this models, are essential for getting the right 
behaviour of the form factors at high invariant 
masses of the meson pairs. 

The ansatz for the pion form factor reads 



OO 



fit 



_n— 6 



(1) 



J dQCD 



with the parameters of the first six radial excita- 
tions fitted to the data and others taken as as- 
sumed in the original model [23] . For details we 
refer the reader to [17] . In the original model the 
coupling constants c^^^ are real, but it is impossi- 
ble to fit the data using real couplings, even if one 
allows that they are different from their model 
values. Thus they were allowed to be complex, 
keeping the normalisation (sum of all couplings 
— 1) unchanged. With these assumptions the fit 
is very good x^/d-o.f = 271/270. This might 
be an indication that the radial p excitations de- 
cay to the same final states, as in this case the 
mixing between them can generate complex cou- 
plings after diagonalisation of the mass matrix. 
This subject was not studied in details as the ex- 
perimental knowledge about the decay modes of 



the p mesons is very poor with the exception of 

p(77o) m- 

The ansatz for the kaon form factors read 



-Ffc+,o(s) — ax+o 



No 



fit 



n=N„ + l 



dQCD 



.n=0 



fit 



dQCD 



fit 



dQCD 



with 



0-K+ — 2 ' "A'O 



(2) 



(3) 



Two versions of the model were used, one with 
an additional assumption Cp^ = Ct^„ (the 'con- 
strained' model) and one without this assump- 
tion (the 'unconstrained' model). TV^, Nu:: Np 
indicate the highest radial resonance used in the 
fit. For the resonances above these values their 
model parameters were used. N^j, = 2, Np — 3 
for the constrained model and = 1, Np = 2, 
N^ = 2 for the unconstrained model. The uncon- 
strained model fits the data better [17], but the 
quality of the second fit is also reasonable. It is 
interesting to observe that the models give dis- 
tinct predictions in the region were there are no 
data available (Fig. [T]), even if in the region were 
data are available their behaviour is very similar. 

It is also of interest to see how much impor- 
tant is the infinite tower of the resonances in the 
form factors. In Fig. [2] the relative difference 
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Figure 1. The experimental data I25I26I27I28I29] 

compared to the model fits results (see text for 
details). 



Figure 2. Comparison of the kaon form factors 
with and without ('no tail') infinite tower of reso- 
nances, c. denotes the constrained model, u. the 
unconstrained one. 



of the form factors without ('no tail') and with 
the infinite tower of resonances is shown. In the 
latter case only Ni ^ \ , i — ^,>xi,p resonance 
contributions are included. As one could guess, 
the contributions of the higher resonances to the 
form factors are more pronounced for higher in- 
variant masses, but they are not negligible even 
in the lower invariant mass range, where the fit- 
ted resonances are situated. A similar behaviour 
of the pion form factor was observed in 1171 . 

3. Description of narrovir resonances 

The narrow resonance (only J/'0, V'(2S') and <\) 
were considered) contributions were included |17| 
into the Monte Carlo generator PHOKHARA for 
the muon, the charged pion and the kaon (both 
charged and neutral) pair production through the 
following substitution Y-ka{s) ~^ ^i?.p(^) ' ^ ^ 
,K+,K° , R = J/il),^{2S) , where 



//,7r. 



r^vp 



is) 



3Ti 



1 — Aa{s) am^ 



Syi Ff (1 + 4) 



s - Ml 



(4) 



(5) 



For P — n and P = tt, Cp = (it was assumed 
that the narrow resonances do not decay directly 
into fJ-'^ fJ,~ and n~^TT~). The (j) contributions to 
the kaon pair production are included in the kaon 
form factor, hence Sk — 0, while Sp — 1 for P — 
/i and P — TT. The notation and the detailed 
description of the narrow resonance contribution 
to the amplitude can be found in [19] together 
with numerical values of the Cp couplings. 

4. PHOKHARA 7.0 

The new form factors described in Section 
[5] and the couplings of tt+tt", K^K~ 

and K'^K^ pairs to narrow resonances described 
shortly in Section [3] were implemented into a 
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Monte Carlo event generator PHOKHARA. It's 
new release (PHOKHARA 7.0) contains also the 
implementation of the 4-pion hadronic currents 
as described in (16) . not included in the pre- 
vious releases of the code. The code can be 



found on the web page |http://if ic.uv.es/~ 
rodrigo/phokhara/ , where also a user guide is 
provided. 

The variance reduction through standard 
change of variables using the Breit-Wigner shape 
of resonances was enough to allow for a reason- 
able acceptance rate. Other technical problems 
with the treatment of the narrow resonances and 
their solution were discussed at length in [17 5^ 
and will not be repeated here. A typical invari- 
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Figure 3. A comparison of the tt+tt^ 'true' (solid 
line) invariant mass spectra and the tt+tt" invari- 
ant mass spectra seen in a detector with a resolu- 
tion of 14.5 MeV ( as for the BaBar detector [3T]) 
- dotted line. 



ant mass spectrum of the final meson or muon 
pair in the vicinity of a narrow resonance, as seen 
in a radiative return experiment, resembles more 
a resolution curve of a detector then the shape 
of the resonance. A typical picture is shown in 
Fig. [21 where the generated mass spectrum of 
the TT+TT^ pair was convoluted with a Gaussian 



with a detector resolution of the BaBar detector 
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Figure 4. The K^K^ invariant mass spectra 
seen in a detector with a resolution of 14.5 MeV 
( as for BaBar detector [3T]) including only ISR 
NLO corrections (ISRNLO) and adding also FSR 
NLO corrections (FSRNLO). 



For the V'(2S') resonance the smearing effects to 
large extend dilute the influence of the resonance 
on the invariant mass spectrum, as shown in Fig. 
|4] for the charged kaon pair invariant mass distri- 
bution. Thus to study it with the radiative return 
method fairly large statistics is required and it is 
difficult to expect that the method can compete 
with the accuracy of scan experiments. 

The contributions from the J/'0 to all stud- 
ied invariant mass spectra is sizable and one can 
study the J/il^ decays with a good accuracy. Thus 
a legitimate question is how big are the contri- 
butions coming from the final state photon ra- 
diation (FSR). If the detector smearing effects 
are not taken into account, the distortions of 
the meson (muon) pair invariant mass distribu- 
tion by the FSR effects, as compared to ISR, 
are huge as shown in Figs. [5] and [51 They 
are still sizable when one takes into account the 
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Figure 5. Comparisons of the muon pair in- 
variant mass distributions taking into account 
only ISRNLO contributions and the complete 
(ISR+FSR)NLO resuh. 



Figure 7. The relative ratio of the muon pair 
invariant mass distributions taking into account 
only ISRNLO contributions and the complete 
(ISR+FSR)NLO result. Detector smearing ef- 
fects are taken into account. 
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Figure 6. Comparisons of the charged kaon 
pair invariant mass distributions taking into ac- 
count only ISRNLO contributions and the com- 
plete (ISR-HFSR)NLO result. 



smearing effects, as evident from Figs. [7] and 
[51 Moreover they do depend on the event selec- 
tion [30] . However in the integrated cross section 
to large extend the corrections cancel. For the 
muon pair production (Tisrnlo = 6.901 pb and 
ciFSBNLO — 6.954 pb, while for the charged kaon 
pair production oisrnlo = 2.450 • 10^^ nb and 
criFSNLO = 2.442 • 10~^ nb. Thus when one uses 
only the integrated spectra the FSR corrections 
can be safely neglected. 

5. Summary 

The implementation into the PHOKHARA 
Monte Carlo event generator of the narrow res- 
onance contributions is presented and the possi- 
bility of studies of the narrow resonances at B 
factories is discussed. 
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